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ABSTRACT 

Context. Our present understanding of high-mass star formation still remains very schematic. In particular, it is not yet clear how 
much of the difference between low-mass and high-mass star formation occurs during the earliest star formation phases. 

Aims. The chemical characteristics of massive cold clumps, and the comparison with those of their low-mass counterparts, could 
provide crucial clues about the exact role that chemistry plays in differentiating the early phases of low-mass and high-mass star 
formation. Water, in particular, is a unique probe of physical and chemical conditions in star-forming regions. 

Methods. Using the HIFI instrument of Herschel we have observed the ortho-NH 3 (1 0 - 0 0 ) (572 GHz), ortho-H 3 0 (1 10 - loi) 
(557 GHz) and N 3 H + (6 - 5) (559 GHz) lines toward a sample of high-mass starless and proto-stellar clumps selected from the 
“Herschel Infrared Galactic Plane Survey” (Hi-GAL). We compare our results to previous studies of low-mass and high-mass proto- 
stellar objects. 

Results. At least one of the three molecular lines was detected in 4 (out of 35) and 7 (out of 17) objects in the £ = 59° and 
(. = 30° galactic regions, respectively. All detected sources are proto- stellar. The water spectra are complex and consist of several 
kinematic components, identified through a Gaussian decomposition, and in a few sources inverse and regular P-Cygni profiles have 
been detected. All water line profiles of the £ = 59° region are dominated by a broad Gaussian emission feature, indicating that the 
bulk of the water emission arises in outflows. No such broad emission is detected toward the £ = 30° objects. The ammonia line in 
some cases also shows line wings and an inverse P-Cygni profile, thus confirming that NH 3 rotational transitions can be used to probe 
the dynamics of high-mass star forming regions. Both bolometric and water line luminosity increase with the continuum temperature. 
Conclusions. The higher water abundance toward the £ = 59° sources, characterized by the presence of outflows and shocks, supports 
a scenario where the abundance of this molecule is linked to the shocked gas. Various indicators suggest that the detected sources 
toward the £ = 30° region are in a somewhat later evolutionary phase compared to the £ - 59° field, although a firm conclusion is 
limited by the small number of observed sources. We find many similarities with studies carried out toward low-mass proto-stellar 
objects, but there are indications that the level of infall and turbulence in the high-mass protostars studied here could be significantly 
higher. 
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1. Introduction 

High-mass stars (OB spectral type, M <; 8 M 0 ), though much 
rarer than low-mass and intermediate-mass stars, are fundamen¬ 
tal in the evolution of galaxies. But, while a reasonably ro¬ 
bust evolutionary sequence has been established for the forma¬ 
tion and early evolution of individual, our present understand¬ 
ing of high-mass star formation still remains very schematic 
(Zinnecker & Yorke 2007, Tan et al. 2014). In particular, some 
of the main issues regarding high-mass star formation that still 
need to be studied are: (i) what is the physical and dynamical 
evolutionary path of massive starless 1 clumps 2 toward the proto- 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im¬ 
portant participation from NASA. 

1 As identified through the absence of Mid-IR emission and/or the 
strength of the PACS IQ pm emission. 

2 In the literature, the term “core” or “fragment” usually refers to 

small-scale objects (S= 0.1 pc or even «0.1 pc), possibly corresponding 


stellar and stellar phases? And, (ii) how much of the difference 
between low-mass and high-mass star formation occurs during 
the earliest star formation phases? 

The high-mass analogues of low-mass Class 0 protostars and 
starless clumps (e.g., Motte et al. 2007, Bontemps et al. 2010, 
Olmi et al. 2010, Olmi et al. 2013) constitute a fundamental 
link between the global processes that regulate star formation 
on large scales and the properties of newly formed (proto)stars 
and clusters. Therefore, to better define the evolutionary path 
from molecular clumps to proto-stars it is crucial to determine 
the physical and dynamical (and, possibly, chemical) properties 
of these early phases of high-mass star formation. 

The formation of both low-mass and high-mass stars from 
the initial dense clumps is accompanied by large changes in the 


to a later stage of fragmentation, while the term “clump” is generally 
used for a somewhat larger (<5 0.1 — 1 pc), unresolved object, possibly 
composed of several cores. The term clump is thus more appropriate to 
refer to the compact objects described in this work. 
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Table 1: Observed transitions. 


Species 

Frequency" 

[GHz] 

Transition* 

£1 

[K] 

£ u 

[K] 

log(A ul ) 

[s’ 1 ] 

ortho-NH^ 

572.4982 

O 

O 

1 

O 

II 

0 

27 

-2.80 

ortho-ICO 

556.9359 

O 

1 

O 

II 

34 

61 

-2.46 

N 2 H + 

558.9665 

7 = 6-5 

67 

94 

-1.96 


Notes. <a> Frequencies from http: //www. cv.nrao. edu/php/splat/ <b> Critical densities are n a ^ 10 7 - 10 s cm 3 for all three lines. 


temperature and density structure in their surrounding envelopes. 
These changes in turn open up different chemical routes. Current 
observational evidence indicates that significant amounts of 
volatile molecules, for instance CO, H 2 O and NH 3 , reside on the 
surfaces of dust grains in dense clumps. These surface molecules 
can be formed either in the gas phase via ion-molecule reac¬ 
tions followed by accretion onto the dust particles, or via direct 
production on the grain surfaces. Different physical conditions, 
for example during the early formation phases of low-mass and 
high-mass stars, could thus determine which of these different 
chemical routes to molecule formation are active and/or domi¬ 
nant. 
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Table 2: Selected observational and physical properties of the source sample: lcr noise level, SSB continuum intensity (T c ), beam¬ 
filling factor (/7bf), evolutionary flag, molecular hydrogen volume density (nffT)), and column density (AYHt)). Luminosity (L), 
mass (M), dust temperature 7’<|, source diameter (Full Width at Half Maximum, or FWHM) and distance (cl) are all derived from the 
Hi-GAL observations. Sources with their name in boldface have at least one detection and thus they could be assigned a Vi sr . 


Source 

Gal. Long, 
[deg] 

Gal. Lat 
[deg] 

L 

[L 0 ] 

M 

[M 0 ] 

T& 

[K] 

Source FWHM 
["] 

d 

[kpc] 

r T' a 

[mK] 

1 CT b 

[mK] 

fi sr C 

[km s _1 ] 

lbf d 

«(H 2 ) e 
[xlO 4 cm“ 3 ] 

mn 2 y 

[xlO 22 cnr 2 ] 

Evol.Flagf 

159-251 

57.984803 

0.0927620 

27.2 

165.0 

10.1 

17.3 

7.0 

0 

9 


0.2 

2.7 

4.82 

1 

159-139 

58.048693 

0.1889100 

117.5 

2118.0 

9.2 

27.8 

9.1 

0 

8 


0.4 

3.7 

14.0 

1 

159-376 

58.171557 

0.7804730 

7.7 

118.0 

8.9 

25.3 

2.6 

0 

8 


0.3 

12.0 

12.0 

1 

159-199 

58.256292 

0.1431680 

67.5 

323.0 

11.5 

31.5 

7.2 

0 

8 


0.4 

0.81 

2.7 

1 

159-203 

58.464407 

0.2996740 

11.5 

81.0 

10.8 

20.0 

3.5 

0 

8 


0.2 

6.5 

6.9 

1 

159-129 

58.464903 

0.4929190 



12.6 

17.9 


0 

8 


0.2 



1 

159-185 

58.512225 

0.5060260 

5.2 

132.0 

8.3 

25.3 

2.8 

0 

8 


0.3 

11.0 

11.0 

1 

159-102 

58.517332 

0.4718170 

8.7 

100.0 

9.2 

25.3 

2.8 

0 

9 


0.3 



1 

159-78 

58.558761 

0.3525380 

33.7 

3296.0 

7.0 

25.3 

6.2 

0 

9 


0.3 

25.0 

58.0 

1 

159-314 

58.569252 

0.0207140 

27.6 

100.0 

11.1 

24.8 

5.7 

0 

10 


0.3 

0.53 

1.1 

1 

159-391 

58.606871 

0.7818820 

25.0 

109.0 

10.7 

25.3 

4.7 

0 

9 


0.4 



1 

159-80 

58.646263 

0.5462850 

11A 

1027.0 

9.1 

25.3 

6.1 

24 

10 


0.3 

8.0 

19.0 

1 

159-345 

58.693246 

0.7685080 

25.7 

149.0 

10.3 

25.3 

4.9 

0 

8 


0.3 

4.5 

10.0 

1 

159-136 

58.754808 

0.6835960 

23.1 

446.0 

8.6 

25.3 

5.5 

0 

8 


0.3 

4.8 

9.9 

1 

159-231 

58.755031 

0.7204550 

88.4 

137.0 

13.8 

33.8 

5.6 

0 

9 


0.5 

0.58 

1.7 

1 

159-293 

58.810992 

0.6553700 

14.5 

3.0 

19.6 

17.8 

5.2 

0 

9 


0.2 

0.11 

0.15 

0 

159-667 

58.865456 

-0.212801 

653.6 

799.0 

12.5 

17.0 

13.7 

0 

8 


0.2 

1.8 

6.3 

2 

159-230 

59.051623 

-0.229814 

23.0 

79.0 

12.1 

27.0 

2.6 

0 

8 


0.4 

6.8 

7.1 

1 

159-57 

59.151275 

-0.268629 

185.7 

281.0 

13.9 

24.2 

6.2 

22 

9 


0.3 

2.4 

5.4 

1 

159-191 

59.208976 

0.1322060 

62.8 

239.0 

11.0 

25.3 

6.7 

0 

8 


0.4 

1.1 

1.0 

1 

159-263 

59.212680 

0.1586920 

14.0 

170.0 

9.2 

25.3 

4.0 

0 

9 


0.5 

0.66 

2.1 

1 

159-150 

59.300967 

0.2429840 

47.3 

210.0 

10.7 

25.3 

5.9 

0 

8 


0.2 

13.0 

20.0 

1 

159-104 

59.304277 

0.2083220 

17.0 

2.0 

23.3 

25.3 

2.7 

0 

9 


0.3 

0.18 

0.18 

1 

159-371 

59.309467 

-0.397915 

15.8 

227.0 

9.6 

19.3 

6.2 

0 

8 


0.2 

3.8 

6.8 

1 

159-31 

59.541533 

-0.101984 

20.7 

5.0 

19.4 

25.3 

2.3 

0 

9 


0.2 

60.0 

30.0 

1 

159-441 

59.633104 

-0.191683 

2259.6 

294.0 

21.6 

15.3 

2.2 

110 

8 

27.8 

0.1 

230.0 

110.0 

2 

159-445 

59.635570 

-0.186242 

787.9 

82.0 

20.9 

23.0 

2.3 

75 

8 

27.4 

0.3 

16.0 

13.0 

2 

159-282 

59.707751 

0.5772220 

30.1 

112.0 

10.9 

25.3 

4.3 

0 

8 


0.1 



1 

159-409 

59.737176 

0.5263750 

60.7 

123.0 

13.2 

27.8 

6.2 

0 

8 


0.4 

0.7 

1.8 

1 

159-444 

59.789451 

0.6307310 



21.8 

9.0 


37 

8 

34.3 

0.1 



2 

159-442 

59.832526 

0.6718730 



25.3 

13.9 


28 

8 

34.2 

0.1 



2 

159-359 

59.862726 

-0.210802 

55.3 

429.0 

9.8 

25.3 

6.4 

0 

11 


0.5 

5.1 

18.0 

1 

159-101 

59.880982 

-0.102392 

85.9 

2814.0 

8.4 

25.9 

6.6 

0 

8 


0.3 

16.0 

41.0 

1 

159-69 

59.891888 

-0.112074 

27.8 

92.0 

12.2 

36.2 

1.9 

0 

8 


0.5 

7.5 

8.0 

1 

159-198 

59.946964 

0.0610780 

5.9 

5.0 

15.3 

13.5 

1.9 

32 

8 


0.1 

8.3 

3.2 

1 


Notes. (a) SSB continuum intensity at 557 GHz. Rms noise at a resolution of 1 km s' 1 with subtracted baseline . (c) Velocity is specified only for 
sources with at least a detection. (rf) Beam-filling factor, see Section 3.4. w See Section 3.4 for a description of these parameters . i ft 0: no detection 
at 70^m, Larson’s criterion (Larson 1981) not fulfilled (starless, unbound). 1: no detection at 70^m, Larson’s criterium fulfilled (pre-stellar). 2: 
detection at 70 fim (proto-stellar). 
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Table 2: continued 




Source 

Gal. Long, 
[deg] 

Gal. Lat 
[deg] 

L 

[L 0 ] 

M 

[Mo] 

Td 

[K] 

Source FWHM 
["] 

d 

[kpc] 

t a 

[mK] 

1 (T b 

[mK] 

v, sr c 

[km s -1 ] 

rib/ 

«(H 2 ) e 

[x 10 4 cm -3 ] 

iV(H 2 T 
[xlO 22 cm' 2 ] 

Evol.Flag^ 

130-376 

29.002507 

0.0689530 

9328.5 

1072.0 

19.8 

21.8 

8.9 

23 

8 

70.5 

0.3 

4.1 

12.0 

2 

130-512 

29.235020 

-0.046735 

1312.0 

1860.0 

13.8 

12.6 

9.9 

63 

8 

61.4 

0.1 

27.0 

51.0 

2 

130-313 

29.862299 

-0.043707 

53790.2 

1866.0 

22.9 

12.2 

8.5 

68 

8 

101.5 

0.1 

47.0 

73.0 

2 

130-110 

29.868288 

-0.007673 

313.8 

1311.0 

10.6 

14.7 

8.5 

27 

8 


0.3 

10.0 

20.0 

1 

130-327 

30.009612 

-0.273200 

12371.5 

964.0 

23.5 

9.2 

6.7 

98 

8 

102.6 

0.1 

120.0 

110 

2 

130-281 

30.121077 

-0.548844 

435.5 

628.0 

12.8 

25.3 

8.1 

44 

8 


0.3 

2.1 

6.4 

1 

130-375 

30.197099 

0.3101250 

15852.5 

2921.0 

17.7 

12.7 

14.2 

0 

9 


0.1 

14.0 

39.0 

2 

130-432 

30.447046 

-0.024628 

9848.3 

271.0 

26.7 

28.0 

11.6 

44 

8 


0.2 

6.5 

22.0 

2 

130-103 

30.575773 

0.1188560 

186.4 

878.0 

10.6 

25.3 

5.7 

44 

8 

55.6 

0.3 

8.4 

18.0 

1 

130-403 

30.627182 

-0.063627 

1739.5 

711.0 

15.8 

28.0 

5.6 

38 

8 


0.4 

5.3 

12.0 

2 

130-73 

30.658809 

0.0450200 

143.7 

962.0 

11.1 

17.7 

5.6 

34 

8 


0.2 

28.0 

42.0 

1 

130-111 

30.661139 

0.1419250 

392.9 

249.0 

14.8 

25.3 

5.6 

35 

10 


0.4 

7.9 

20.0 

1 

130-132 

30.675761 

-0.166004 

1223.9 

1559.0 

14.3 

27.4 

9.0 

0 

8 


0.4 

3.0 

11.0 

1 

130-42 

30.693536 

-0.148593 

1082.2 

723.0 

16.7 

17.4 

5.7 

50 

9 

90.4 

0.2 

22.0 

32.0 

2 

130-304 

30.703241 

-0.068204 

49112.4 

3612.0 

22.5 

20.4 

6.0 

350 

12 

90.7 

0.2 

57.0 

100.0 

2 

130 - 43 ? 

30.822869 

0.1343020 

2405.6 

16437.0 

10.2 

32.6 

8.8 

50 

8 

108.4 

0.5 

20.0 

86.0 

2 

130-230 

30.846821 

0.3426510 

186.5 

425.0 

12.0 

25.3 

6.2 

0 

8 


0.4 

5.7 

21.0 

1 


Notes. (g> Tentative detection (see Section 3.2). 
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The chemical characteristics of massive cold clumps, and the 
comparison with those of their low-mass counterparts, are only 
starting to emerge (Russeil et al. 2010, Marseille et al. 2010). 
While governed by the same processes as lower mass objects, 
the much larger luminosity of massive stars during all their evo¬ 
lutionary phases will greatly affect the outcome. Therefore, the 
comparison of chemical models (e.g., with different weights of 
grain vs. gas-chemistry) with aimed spectral line observations 
will provide crucial clues about the exact role that chemistry 
plays in differentiating the early phases of high-mass star for¬ 
mation from those of low-mass star formation. 

The (1 q - Oo) ground-state rotational transition of 0 -NH 3 
(e.g.. Ho & Townes 1983) has an upper state energy of 27 K 
(see Table 1), comparable to the rotation-inversion transitions at 
cm wavelengths. The critical density is, however, four orders of 
magnitude higher and is thus a better probe of the physical con¬ 
ditions of the dense high-mass starless clumps detected by the 
Hi-GAL survey. Moreover, its ortho symmetry form makes it 
a complementary tool to the commonly used rotation-inversion 
transitions of para-NH 3 observable from the ground. Both NH 3 
and N 2 H + are well known tracers of low-mass pre- and proto- 
stellar clumps. 

Water is another important component in the chemistry of 
regions of star formation. In cold gas, water is mainly produced 
on dust grains, and released to the gas phase when the temper¬ 
ature increases (^ 100 K). Water is thus not expected to to be 
found with high abundances in cold, starless clumps (e.g., Bergin 
& Tafalla 2007, Caselli et al. 2012, Wirstrom et al. 2014). In 
fact, recent results from the Herschel Key Program “Water in 
star-forming regions with Herschel” (WISH; van Dishoeck et al. 
2011 ) have shown that during the quiescent phases of star for¬ 
mation, the water abundance is very low. The embedded stage 
of star formation, however, shows very broad and complex line 
profiles consisting of many dynamical components (Kristensen 
& van Dishoeck 2011, Kristensen et al. 2012). The WISH team 
also found, comparing typical low-, intermediate-, and high- 
mass young stellar objects, that their spectra were remarkably 
similar, characterized by broad, complex line profiles and con¬ 
sisting of multiple components. Our observations mainly con¬ 
sists of high-mass clumps in various evolutionary stages, and 
thus we will also try to identify possible evolutionary differ¬ 
ences. 

The outline of the paper is as follows. In Section 2 we dis¬ 
cuss how our sources were selected and observed with Herschel. 
In Section 3 we describe how the main kinematical and phys¬ 
ical parameters of the detected sources were derived. The gas 
kinematics is further discussed in Section 4 and we draw our 
conclusions in Section 5. 

2. Observations and data reduction 

2.1. Source selection 

We have selected our targets from the preliminary source cat¬ 
alog obtained toward the two Galactic fields, at £ = 30° and 
£ = 59°, observed for the “Herschel Infrared Galactic Plane 
Survey” (Hi-GAL, Molinari et al. 2010) during the Science 
Demonstration Phase of Herschel (Elia et al. 2010). The Hi- 
GAL (and Hi-GAL360) Project used SPIRE and PACS in paral¬ 
lel mode to carry out an unbiased imaging survey of the Galactic 
Plane, uniformly sampling a 2°-wide strip, with -60° < t < 
60°, in the 70, 170, 250, 350 and 500/rm photometric bands. 
The source-identification, flux-extraction and spectral energy 
density-building in these two Hi-GAL maps were carried out by 


adapting to the Herschel SPIRE/PACS wavebands the methods 
and techniques described by Netterfield et al. (2009) and Olmi 
et al. (2009). From the preliminary catalogs of the ( = 30° and 
£ = 59° fields, we have extracted both high-mass starless clumps 
and high-mass proto-stellar objects. 

Candidate high-mass starless clumps have been selected 
by imposing that M c i ump > 100 M Q (M c i ump > 200 M 0 ) and 
Lfir < 100 L q (Lfl r < 500 L q ) for the £ - 59° {£ = 30°) field. 
Moreover, we have checked that no compact MIPS (“Multiband 
Imaging Photometer for Spitzer”, Rieke et al. 2004) 24 /un 1 emis¬ 
sion were present at the position of the candidate high-mass star¬ 
less clump. Candidate high-mass proto-stellar objects do have 
instead a compact MIPS 24 /urn counterpart. Finally, the sources 
have been visually inspected to reject ambiguous detections. 
This procedure resulted in a total of 81 clumps, including the 
four candidate high-mass starless clumps already observed with 
the VLA and the Effelsberg 100-m telescope in the NH 3 G, 1) 
and (2,2) inversion lines by Olmi et al. (2010). Out of the initial 
sample, 52 sources were finally observed with Herschel'. 17 ob¬ 
jects in the £ = 30° region and 35 in the £ = 59° region. The 
observed sources and their basic physical parameters are listed 
in Table 2. Sources with a detection in at least one of the three 
observed lines have been highlighted in boldface. 

2.2. Herschel observations 

Using the HIFI instrument (de Graauw et al. 2010) of Herschel 
(Pilbratt et al. 2010) we have observed the 0 -NH 3 G 0 - Oo) 
line at 572 GHz in the upper side band of band lb, and simul¬ 
taneously in the lower side band both the 0 -H 2 OG 10 - loi) 
line at 557 GHz, and the N 2 H + (6 - 5) line at 559 GHz (see 
Table 1). We used the wide band spectrometer with a bandwidth 
of 4 GHz and an effective spectral resolution of 1.1 MHz (AV = 
0.27 km s 1 ). The half-powerbeam width and the main-beam ef¬ 
ficiency at 572 GHz are 36 "and 0.62, respectively (Mueller et 
al., 2014 3 ). The total calibration uncertainties are - 9 % for band 
l 4 (Roelfsema et al. 2012). We have used the dual beam switch¬ 
ing (DBS) observing mode. The reference beams in these obser¬ 
vations were located within 3' on either side of the source. The 
observations were performed during the Herschel Open Time 
(OT2) on April 19 and 20, 2013, and the corresponding observa¬ 
tion identifier (OBSID) range is [1342270449,1342270604], 

The data were processed using the standard Herschel 
Interactive Processing Environment 5 , version 11.1, up to level 2, 
providing fully calibrated double side band spectra in the T* K an¬ 
tenna temperature intensity scale where the lines are calibrated 
on a single side band (SSB) scale. Because HIFI is intrinsically a 
double side band instrument, the continuum has to be divided by 
two to be properly scaled. The data quality is excellent with very 
low intensity ripples in most cases, typically below a few per¬ 
cent of the double side band continuum. However, in some cases 
the baselines are bad enough that the large rms (<: 20 - 30 mK) 
does not allow a reliable estimate of the continuum temperature 
to be obtained (see Section 3.3). Therefore, for these sources 
we set the continuum temperature to 0. The FITS files were ex¬ 
ported to the spectral line analysis software packages xs 6 and 


3 http://herschel.esac.esa.int/twiki/pub/Public/ 
/HifiCalibrationWeb/HifiBeamReleaseNote_Sep2Q14.pdf 

4 http://herschel.esac.esa.int/Docs/HIFI/html/ch5.html 

5 http://herschel.esac.esa.int/HIPE_download.shtml 

6 http://www.chalmers.se/rss/oso-en/observations/ 
data-reduction-software 
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Fig. 1: Top panel. Gaussian decomposition of the o-H 2 0(1io - 
lot) spectrum towards the source 159-441. The individual 
Gaussian components are shown as solid color lines overlaid on 
the (continuum-subtracted) spectrum. The overall Gaussian fit is 
shown as a red dashed line, and the residual deviation at a ve¬ 
locity of - 55 km s-' is visible. The intensity of the component 
in absorption has been multiplied by a 0.2 factor. Bottom panel. 
Overlay of the 0 -H 1 O (lio-loi), 0 -NH 3 (l 0 -0 0 ) andN 2 H + (6-5) 
lines towards the source 159-441. 


CLASS 7 . Both polarisations and all three LO-settings were in¬ 
cluded in the averaged noise-weighted spectra. The resulting av¬ 
erages were convolved to a resolution of 1 km s ’ 1 in sources 
with a low signal-to-noise ratio. All HIFI spectra, in T* A units. 


are shown in Appendices A (C = 30° region) and B (£ — 59° 
region). 


3. Results 

3.1. Distance determination 

Assigning a distance to sources detected with a photometer is a 
crucial step in giving physical significance to all information ex¬ 
tracted from the Hi-GAL data. While reliable distance estimates 
are available for a limited number of known objects (e.g., Hii re¬ 
gions, see Russeil 2003, and masers, see e.g. Green & McClure- 
Griffiths 2011), this information does not exist for the majority 
of Hi-GAL sources. We therefore adopted the scheme presented 
by Russeil et al. (2011) to assign kinematic distances to large 
lists of sources. A 12 CO (or 13 CO) spectrum (e.g., from the BU- 
FCRAO Galactic Ring Survey, or GRS, Jackson et al. 2006) is 
extracted in the line of sight of every Hi-GAL source. Then, the 
velocity, Vj sr , of the brightest spectral component is assigned to 
that specific source thereby allowing calculation of a kinemati- 
cal distance. By using extinction maps (derived from the 2MASS 
point source catalog, see e.g. Schneider et al. 2011) and a cata¬ 
log of sources with known distances, such as Hit regions, masers, 
and others, the distance ambiguity is resolved and a recommen¬ 
dation given. In this way it is possible to produce a “distance 
map” having the same pixel size of the CO cube as used to ex¬ 
tract the spectra for every target, where the value of the pixel 
is the assigned distance of the Hi-GAL source(s) falling in that 
pixel. The typical error on the kinematic distance 8 (column 8 in 
Table 2) is estimated to vary between ~ 0.6-0.9 kpc in the range 
of longitude ( ~ 30° - 60°. 

3.2. Profile decomposition 

The observed line profiles, and in particular those of the water 
spectral lines, show complex shapes, and likely represent contri¬ 
butions from proto-stellar envelopes, molecular outflows, fore¬ 
ground clouds as well as infall motions. These contributions, 
both in emission and in absorption, are disentangled and param¬ 
eterized by fitting multiple Gaussians to the line profiles. We 
show an example in Fig. 1 and list all the results in Table 3. 
The Gaussian fits have been performed using up to four velocity 
components and all line parameters were allowed to vary. A sim¬ 
ilar decomposition of the H 2 0(1io~ lot) (557 GHz) line has also 
been used by Kristensen et al. (2012) and Herpin et al. (2012). 

Table 3 shows that the H 2 0 line profiles of the ( = 59° 
sources are very similar. Almost all water spectra can be de¬ 
scribed as the sum of a narrow (FWHM~ 3-5 km s’ 1 ), a medium 
(FWHM~ 6 - 10kms’ 1 ) and a broad (up to - 30kms’ 1 ) ve¬ 
locity component. The narrow component is mostly seen in ab¬ 
sorption, whereas the medium and broad components have pos¬ 
itive intensities. Source 159-444 shows an inverse P-Cygni pro¬ 
file overlaid on a broad component. In source 159-441 and, to 
a lesser extent, in 159-445 (please note that these two sources 
are only 21.5" apart) if three Gaussian components are used in 
the fit the red-shifted wing of the spectrum still shows a resid¬ 
ual deviation Gt 3 <x, see Fig. 1) from the overall fit at a ve¬ 
locity of - 55 km s’ 1 . This residual can be eliminated with an 
additional Gaussian component, which may either have a high 
velocity (<; 50 km s’ 1 ), or a velocity comparable to the other 
components but with a wide line width (<; 4010ns’ 1 ). In the 


7 CLASS is part of the GILDAS software package developed by 8 This estimate assumes that the distance ambiguity has been cor- 

IRAM. rectly solved. 
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Table 3: Gaussian fits of the Herschel-HIFI observations. 


Source 

HiO Comp." 

E/A 6 

H 2 0 

V, s r 

[km s' 1 ] 

n 

[K] 

AV 

[km s' 1 ] 

Vlsr 

[km s' 1 ] 

N,H + 

n 

[K] 

AV 

[km s' 1 ] 

Vlsr 

[km s' 1 ] 

NH, 

n 

[K] 

AV 

[km s' 1 ] 

159-441 

M 

E 

27.9 

0.89 

6.3 

27.8 

0.83 

3.1 

27.7 

0.81 

3.6 


B 

E 

27.0 

0.37 

26.0 




28.2 

0.09 

10.1 


N 

A 

26.9 

-1.1 

3.7 




27.7 

-0.46 

2.0 

159-442 

M 

E 

36.7 

0.24 

6.3 

34.2 

0.14 

3.6 

34.8 

0.13 

6.4 


B 

E 

36.7 

0.12 

23.6 








N 

A 

34.9 

-0.30 

4.8 




35.1 

-0.05 

2.5 

159-444 

N 

E 

32.0 

0.21 

2.9 

34.3 

0.20 

4.6 

34.4 

0.21 

5.8 


B 

E 

31.2 

0.08 

20.7 








N 

A 

27.7 

-0.08 

1.0 








N 

A 

33.9 

-0.14 

4.8 




35.4 

-0.18 

3.4 

159-445 

M 

E 

28.0 

0.67 

7.2 

27.4 

0.67 

3.1 

27.4 

0.64 

4.2 


B 

E 

28.4 

0.24 

28.2 








N 

A 

27.0 

-0.80 

3.7 




27.7 

-0.25 

1.7 

130-42 






90.4 

0.03 

1.9 

91.3 

0.08 

4.2 

130-43 rf 

N 

E 

108.4 

0.07 

2.7 


< 0.02" 



< 0.02" 


130-304 






90.7 

0.86 

4.1 

87.8 

0.33 

3.1 


M 

A 

92.9 

-0.35 

8.5 




93.7 

-0.31 

5.7 

130-313 

M 

E 

101.2 

0.07 

6.5 

101.5 

0.08 

2.7 

101.5 

0.14 

5.8 

130-327 






102.6 

0.07 

3.6 





M 

A 

103.8 

-0.11 

5.7 




104.6 

-0.08 

4.6 

130-376 

N 

E 

70.5 

0.07 

3.2 


< 0.02" 


70.5 

0.05 

2.7 

130-512 

M 

E 

63.2 

0.11 

5.5 

61.4 

0.06 

2.9 

61.8 

0.15 

2.5 


N 

A 

59.5 

-0.06 

2.4 








N 

A 

64.0 

-0.06 

1.6 








Notes. ) o) Definition of the Gaussian components of the water line: N=narrow, M=medium, B=broad. <w E=emission, A=absorption. 3a limit. 
Id> Tentative detection (see Section 3.2). 


first case the additional velocity component could be interpreted 
as an extremely high velocity (EHV) component. Similar EHV 
components have previously been reported, e.g., by Kristensen 
et al. (2012) toward several low-mass protostars and by Leurini 
etal. (2014) toward the massive protostar IRAS 17233-3606, and 
are usually associated with shocked gas. However, we do not 
have further evidence (such as spectra of other water lines and/or 
spectra at different positions) to confirm that an EHV component 
indeed exists in source 159-441. In the second scenario, an addi¬ 
tional broad velocity component must be invoked, which may 
be caused by the fact that molecular outflows are not Guassian 
shaped. Due to this ambiguity, we choose to show only three 
Gaussian components for this source. Section 4.2 discusses the 
water line profile components in detail. 

The situation is very different in the { = 30° region, with 
no two profiles being identical. We have no clear detection of a 
broad component in the water spectra, and the line profile can 
be fit by fewer Gaussian components. Ammonia also displays 
complex line shapes with mixed emission and absorption, and 
only weak signs of outflows. In contrast to water and ammonia, 
the N 2 H + spectra do not show any line asymmetries and we ob¬ 
serve only emission. In addition, two sources also show a regular 
or inverse P-Cygni profile, either in the water line (130-512) or 
in the NH 3 spectrum (130-327) (see Section 4.2.5). Tables 2, 3 
and 4 mark source 130-43 as a tentative detection because recent 
(June 2015) follow-up observations with the 20-m telescope of 


the Onsala Space Observatory have detected the NH 3 ( 1 , 1 ) line 
toward this source at a velocity of - 95.3 km s' 1 , instead of 
108.4km s ' 1 as shown in Table 3 and Appendix A. Therefore, 
either the water line in source 130-43 is an artefact of some type, 
or both lines are actually real and thus the issue remains as to 
whether or not they belong to the same object. With our current 
data we cannot answer this question and therefore in this work 
we will assume that the water line detection is real. 


3.3. Correlations between dust- and gas-derived physical 
parameters 

In order to better understand the nature of the detected sources 
in the C = 30° and t - 59° regions, as well as to determine 
the origin of the H 2 0 emission, we have performed correlations 
between several parameters. We begin with a comparison which 
includes both dust-derived physical parameters, as determined 
from our previous Hi-GAL observations, and HIFI-derived pa¬ 
rameters, such as the continuum temperature, T c . Then, the total 
H 2 G integrated line intensity is compared to different parame¬ 
ters. 

In Fig. 2 we plot the clump luminosity vs. mass, as de¬ 
rived from our previous Hi-GAL observations. Since both are 
distance-dependent quantities, only sources where a distance 
could be determined are shown. The plot shows that the clumps 
found in the t = 59° region are characterized by lower lumi- 
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Fig. 2: Clump bolometric luminosity vs. mass, as derived from 
the Hi-GAL observations for all sources where a distance could 
be determined. Sources in the £ = 30° and £ — 59° regions are 
represented by the (black) “+” sign and (red) open squares, re¬ 
spectively. Sources with a detection in at least one of the three 
observed spectral lines are shown with thicker (green) symbols. 
The dashed lines are loci at constant temperature. Roughly or¬ 
thogonal to these are loci (dotted lines) of constant 250 pm flux 
density, ranging from 0.1 to 500 Jy, assuming a modified black- 
body spectral energy distribution with p — 1.5 and a fixed dis¬ 
tance of ~ 6kpc, equal to the median of the distances of the 
observed sources. The black dots and the lines joining them rep¬ 
resent the evolutionary tracks of clumps with starting envelope 
mass equal to 120, 470 and 10 3 M 0 (see Section 4.4). 


nosity and mass compared to those in the t = 30° region (see 
the detailed statistical analysis by Olmi et al. 2013 and Olmi 
et al. 2014a). Sources with a detection in at least one of the three 
observed spectral lines are clearly found at higher luminosity 
(L 2 10 3 L 0 ), as also confirmed by Table 5. In Fig. 2 we also plot 
a series of evolutionary tracks which we discuss in Section 4.4. 

In Fig. 3 we plot the dust-derived temperatures (75, column 
6 in Table 2) and luminosities of the clumps in the ( = 30° and 
( = 59° regions vs. T c . The continuum flux density associated 
with T c is likely to have a contribution from both free-free emis¬ 
sion and warm dust thermal emission. We do not have an alterna¬ 
tive method to determine the contribution of free-free emission 
in our sources. Furthermore, because we have only measured 75 
at one frequency, we cannot determine how T c varies with fre¬ 
quency, which could otherwise be used to discriminate between 
dust and free-free emission, van der Tak et al. (2013), however, 
found that the free-free emission in their massive proto-stellar 
clumps is negligible compared to warm dust emission. If this 
is true also for our sources, then 75 should correlate with the 
source bolometric luminosity. This is indeed shown in the top 
panel of Fig. 3, where we also used the Bayesian IDL routine 
LINMIX_ERR to perform a linear regression to find the slope of 
the best fit line. The source luminosity is expected to increase 
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Fig. 3: Top panel. Clump bolometric luminosity vs. continuum 
temperature, 75. The solid line represents the linear fit to all 
points (see text); the Spearman rank coefficient is 0.65. Bottom 
panel. Dust-derived temperature, 75, vs. continuum temperature. 
Symbols are as in Fig. 2. 


with time, as the clump evolves from the starless to the proto- 
stellar phase (see, e.g., Molinari et al. 2008). Therefore, the top 
panel of Fig. 3 is also likely to represent an evolutionary se¬ 
quence. However, the bottom panel of Fig. 3 shows that there is 
only a marginal correlation between 75 and 7' c , which are both 
distance-independent. The large scatter is certainly due in part 
to the large uncertainties in 75, since the baseline rms was in 
general quite large (2 20 - 30 mK). 

In Fig. 4 we plot various line integrated intensities. The in¬ 
tegrated area of the profile has been determined by summing 
all components in emission. For a Gaussian line profile, the in¬ 
tegrated emission is 7'J[ k AV y/n/(4 In 2), where is the peak 
antenna temperature and AF is the FWHM of the line. Sources 
with pure absorption spectra in any of the plotted spectral lines 
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have not been considered. Despite the low number of data points, 
the top and middle panels show an excellent correlation between 
the integrated intensities of the H 2 O and NH 3 lines, and also be¬ 
tween those of the two N-bearing molecules. In the middle panel 
source 130-304 has not been included in the plot because it has a 
completely self-absorbed water line profile, and the NH 3 line has 
an inverse P-Cygni profile. Therefore, the integrated intensity (in 
emission) of ammonia is very uncertain. The reasonable agree¬ 
ment between the integrated intensities of ammonia and N 2 I-D 
suggests that they are emitted from similar volumes of gas, as¬ 
suming similar excitation conditions. 

Finally, in the bottom panel of Fig. 4 we plot the quantity 
f 7 ’^[H 20 ] dv x d 2 . proportional to the H 2 O line luminosity, vs. 
the continuum temperature, T c . Although we do not show any 
linear fit in this case, there seems to be a general trend of in¬ 
creasing line luminosity with T c . The 130 source with the high¬ 
est water line luminosity 60Kkms~' kpc 2 ) is 130-512 which, 
however, has a P-Cygni profile (see Section 4.2.5) and thus the 
actual line luminosity is uncertain. Except for this point, the plot 
shows that in these objects the “typical” water line luminosity 
toward the l = 30° sources is lower compared to that of the 
l = 59° region. 


3.4. Column densities and molecular abundances 


3.4.1. RADEX modelling 

Since we lack direct information about the excitation temper¬ 
atures and opacities of the observed transitions, we use the 
non-equilibrium homogeneous radiative transfer code RADEX 9 
(van der Tak et al. 2007) to correct for possible populations of 
molecules in unobserved excited levels, in order to relate the ob¬ 
served emission and absorption components to the column den¬ 
sities of H 2 0, NH 3 and N 2 H + . 

The principal parameters needed to specify a model are the 
kinetic temperature, 7\, the average number density of molecu¬ 
lar hydrogen «h 2 , and the line width (see for example, Persson 
et al. 2009). The density of each source is computed from the Hi- 
GAL observed size and mass, assuming a homogeneous, spher¬ 
ical cloud. The kinetic temperature is assumed to be equal to the 
dust temperature obtained from our Hi-GAL observations. All 
parameters are listed in Table 2. The background radiation field 
is an average of the Galactic background radiation in the solar 
neighbourhood and of the cosmic microwave background radia¬ 
tion. 

The column density obtained from RADEX is varied until the 
model line brightness temperature, 7b, matches the amplitude 
of the Gaussian fit (Tab. 3). The line brighness temperature is 
obtained from the antenna temperature after correcting for main- 
beam efficiency, rj m b, and beam-filling factor, 77 bf, i.e.: 


e 1 . 

r * rji rrt HID / 1 \ 

A — b 7]mb ^7bf — T b 77mb “2 7>2 (1) 


Lig.4: Top panel. Integrated H20(1io - lot) line intensity plot¬ 
ted vs. the integrated NH 3 (lo - Oo) line intensity. Symbols are 
as in Lig. 2. The solid line represents the linear fit to all points. 
The Spearman rank coefficient is 0.94. Middle panel. Integrated 
NH3(1o - Oo) line intensity vs. the integrated N 2 H + (6 - 5) 
line intensity. Source 130-304 has not been included in the plot 
(see text). The Spearman rank coefficient is 0.96. Bottom panel, 
f r*[H 2 0] dv x d 2 . proportional to the H 2 O line luminosity, vs. 
continuum temperature, T c . Symbols are as in Lig. 2. 


with 6 m b and 9 S being the LWHM of the main-beam and source, 
respectively. The values of r; M are also listed in Table 2. Lor 
the components in absorption, we match the observed and mod¬ 
elled integrated opacity instead of the brightness temperature 
(see Section 4.2.3). 

The computed RADEX results are listed in Table 4. The re¬ 
sults are not very sensitive to changes in density because of 
the high critical densities (n cr ^ 10 7 - 10 8 cm -3 ) of all three 

9 http://www.sron.rug.nl/ vdtak/radex/radex.php 
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Table 4: Column densities and abundances from RADEX modelling 13 . 


Source 

Comp. 6 

Vi sr 

[km s _1 ] 

A(o-H 2 0) c 

[xl0 14 cnr 2 ] 

X(o-H 2 0) rf 
[x 10 9 ] 

v, sr 

[km s _1 ] 

A(N 2 H + ) C 
[xlO 14 cnr 2 ] 

A(N 2 H + )'' 

[xlO- 9 ] 

v, sr 

[km s<] 

A(o-NH 3 ) c 
[ xlO 14 cm<] 

A(o-NH 3 ) j 

[xlO- 9 ] 

159-441 

M 

27.9 

450 

41 

27.8 

0.34 

6.2 

27.7 

37 

3.4 


B 

27.0 

67 





28.2 

1.3 



N 

26.9 

0.34 

0.06 




27.7 

0.05 

0.01 

159-442 

M 

36.7 



34.2 



34.8 




B 

36.7 










N 

34.9 






35.1 



159-444 

N 

32.0 



34.3 



34.4 




B 

31.2 










N 

27.7 










N 

33.9 






35.4 



159-445 

M 

28.0 

350 

270 

27.4 

32 

25 

27.4 

60 

46 


B 

28.4 

300 









N 

27.0 

0.36 

0.6 




27.7 

0.03 

0.04 

130-42 





90.4 

4.3 

1.3 

91.3 

6.0 

1.9 

130-43® 

N 

108.4 

5 

0.58 


<30 

<4 


<0.6 f 

<0.06 f 

130-304 





90.7 

33 

3.3 

87.8 

4.9 

0.49 


M 

92.9 

>2 f 

>0 .¥ 




93.7 

0.52 

0.1 

130-313 

M 

101.2 

15 

2.19 

101.5 

3.1 

0.42 

101.5 

10 

1.4 

130-327 





102.6 

2.4 

0.22 





M 

103.8 

0.73 

0.1 




104.6 

0.21 

0.04 

130-376 

N 

70.5 

39 

33 


<8 

<7 

70.5 

6.1 

5.1 

130-512 

M 

63.2 

110 

22 

61.4 

26 

5.1 

61.8 

30 

5.9 


N 

59.5 

0.3 

0.1 








N 

64.0 

0.2 

0.08 








Notes. (o) Using our derived n( H 2 ) and '1\- = rd ust (Table 2). (w N=narrow component, M=medium, B=broad. Column density corrected for main 
beam efficiency and beam-blling factor (Table 2). Abundance with respect to molecular hydrogen (Table 2). In the absorption components we 
have used N( H 2 )/2 assuming that the absorption is probing approximately half of the total column density in front of the background continuum. 
(c ' 3cr limit. Optically thick; r is assumed to be 5.Tentative detection (see Section 3.2). 


molecular transitions. However, the column density is sensitive 
to changes in kinetic temperature if 7\ $ 15 K. For sources with 
non-detections in one (or two) of the molecular transitions, the 
upper limit to the column density has been estimated by using 
the 3 cr limits for 7’* listed in Table 3. 


Abundances with respect to molecular hydrogen, ATmol) = 
/V(mo])//V(H 2 ), are estimated using the derived volume densities 
and observed source sizes to compute ATH 2 ). For the compo¬ 
nents in absorption we divide the estimated MH 2 ) by a factor 
of two since the absorption is probing half of the total column 
density that is in front of the continuum source. For the outflow 
components we cannot estimate the abundance since ATH 2 ) is 
likely orders of magnitude lower than that of the whole cloud. 


Table 4 indicates that the water molecular abundance is pos¬ 
sibly higher toward sources in the £ = 59° region. Our results 
may be affected by the uncertainties in the derived Hi-GAL 
physical parameters and other assumptions. However, combined 
with the result that the typical water line luminosity toward the 
£ — 30° sources (where no broad component is detected) is lower 
compared to that of the £ — 59° region (see bottom panel of 
Fig. 4), we conclude that our observations support a scenario 
where the water abundance is linked to shocked gas. 


3.4.2. Column density for the absorption component 

We also used an alternative approach to compute the o- H 2 O col¬ 
umn density of the envelope gas responsible for the observed 
absorption in the water lines, which we also use to double-check 
the RADEX results. We first derived the optical depth at the dip of 
the absorption profile as: 


Aiip — — In 



( 2 ) 


where, in order to include the line contribution to the absorp¬ 
tion background, the antenna temperature 7’°'" represents the 
sum of all Gaussian components in emission (at the velocity 
of the dip in the absorption component) and of the continuum 
temperature T c . This procedure assumes that the absorber is lo¬ 
cated in front of both the continuum and the line emitting source. 
7^ lp represents the peak absolute value of the antenna tempera¬ 
ture of the Gaussian component in absorption. The integrated 
optical depth of the absorption conponent is then calculated as 
Tabs = f t dv = 1.06 Tjip AF a bs where AV a b s represents the line 
FWHM of the Gaussian component in absorption. Finally, the 
column density can be estimated as: 


Mot - M = 


87ry ulg' 

c 3 A u \g„ 


1 - exp 


hv iii \ ’ 

kT ex ) 


^"abs 


(3) 
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Table 5: Mean values of bolometric luminosity, mass and continuum temperature. 


Region 

Sources with detection 


Sources with no detection 0 


Mean L 

Mean M 

Mean T c 

Mean L 

Mean M 

Mean T c 


[L e ] 

[M 0 ] 

[K] 

[L 0 ] 

[M e ] 

[K] 

t = 30 

16200 ± 17600 

3430 ± 3300 

93 ±65 

2145 ±2570 

690 ± 305 

37 ±5 

f = 59 

1520 ±735 

190 ± 105 

63 ±30 

90 ±64 

440 ± 390 

26 ±4 


Notes. (a> Only sources where all three parameters could be measured were included. 


where A u i and v u \ represent the Einstein A coefficient and fre¬ 
quency associated with the transition between the lower and up¬ 
per level, indicated by the subscripts / and u, respectively, and 
with degeneracies gi and g„. Physical constants are indicated 
with standard symbols and T ex represents the excitation temper¬ 
ature. Since we do not know T ex we will assume for the mo¬ 
ment that T ex <sc hv u \/k - 26.6 K for the 557 GHz water line. 
Therefore, the exp factor in Eq.(3) can be neglected. 

We thus obtain remarkably similar 0 -H 2 O column densities, 
varying from about 10 13 cm ~ 2 to ^ 5 x 10 l 3 cm ~ 2 (consistent 
with the RADEX results), with the highest value obtained toward 
the source 130-304, where in fact the absorption appears to be 
almost saturated, since the line profile is nearly down to the 
zero-temperature level. The column densities of the outer en¬ 
velope are thus quite similar to those estimated for both low- 
(Kristensen et al. 2012) and high-mass (van der Tak et al. 2013) 
star forming regions. 


4. Discussion 

4. 1. Average physical conditions in the £ = 30° and £ = 59° 
regions 

As shown by Tables 2 and 3, in the £ = 30° region about 41% of 
the sources have a detection in at least one of the observed spec¬ 
tral lines. This percentage drops to about 11% in the C — 59° 
region. We also note that all of the sources with at least one de¬ 
tection are proto -stellar in both galactic regions. The significant 
difference in detections between the two regions may be a conse¬ 
quence of the different fraction of proto-stellar sources observed 
in each region, corresponding to 59% and 14% in £ = 30° and 
£ = 59°, respectively. In fact, the fraction of profo-stellar sources 
with at least one detection corresponds to 70% and 80% in the 
£ — 30° and C = 59° regions, respectively, making the difference 
in the detection rate less significant. 

In Table 5 we list the mean values of the clump mass and 
bolometric luminosity, as derived from the Hi-GAL observa¬ 
tions, for the observed sources in the two galactic regions. We 
note that the mean values of the luminosity for sources with a 
detection in one of the observed spectral lines are much higher 
than the luminosity of sources with no detection at all. The scat¬ 
ter around the mean value is quite high and thus our conclu¬ 
sion is only tentative and will need more data to be confirmed. 
Furthermore, the mean values in the £ = 59° region are derived 
from only two sources, since in sources 159-442 and 159-444 the 
distance could not be determined reliably. 

As with luminosity. Table 5 shows that T c has a higher mean 
value (by more than a factor of 2 ) in sources with at least one de¬ 
tected spectral line, for both observed regions. Given that all of 
these sources are proto-stellar, the mean values of T c for sources 
with and without detection are in fact very similar to the mean 
values of T c for proto-stellar and starless sources, respectively. 
Compared to luminosity and mass, the mean values of T c show a 


lower scatter, which is still quite high only in the case of sources 
with detection in the £ = 30° region. In fact, the scatter for the 
case of sources with no detection is quite low in both regions, 
despite the low number of sources. Given the marginal correla¬ 
tion between luminosity and T c (see Section 3.3) this result is 
not totally unexpected. However, given that the continuum tem¬ 
perature is a distance-independent quantity this result is more 
significant in charaterizing the sources. Finally, from Table 5 we 
also note that the mean values of the mass for sources with and 
without a detection do not show any specific relation. Thus, if we 
assume tentatively that the difference in the mean value of lumi¬ 
nosity (and T c ) between sources with and without detections is 
indeed real, then this might be a consequence of evolutionary ef¬ 
fects rather then an effect associated with more massive clumps 
(see also Section 4.4). 

4.2. Gas kinematics 

Our decomposition of the H 2 O line profiles (see Section 3.2) is 
similar to that of Kristensen et al. (2012) toward low-mass pro¬ 
tostars, and to the description of the H 2 O profiles of van der Tak 
et al. (2013) toward high-mass protostars. Herpin et al. (2012) 
found a similar profile decomposition toward the massive pro¬ 
tostar W43-MM1. In the following we first describe the general 
properties of the line shapes in the £ = 30° and £ = 59° regions, 
and then we discuss the possible physical origin of each compo¬ 
nent. 

4.2.1. Line profile components in the £ = 59° region 

Here we summarize the main observational properties of the line 
profiles of the three transitions, H 20 ( 1 io - lot)- NH 3 G 0 - Oo), 
and N 2 H + (6 - 5) toward the £ = 59° region. 

H 2 O. All detected sources of the £ = 59° region have a 
strikingly similar water line-profile, which is always charater- 
ized by an absorption dip at the source velocity (as given by 
N 2 H + ). This absorption dip does not appear to be saturated. 
As shown by Table 3, the offset between these absorption 
components and the source velocity does not appear to be 
significant. The broad and medium Gaussian components have 
line widths exceeding those of N 2 H + , which is tracing the 
more quiescent ambient gas, and are thus likely originating in 
proto-stellar outflows. The water line profile in source 159-444 is 
somewhat different from that of the other three detected sources 
in the £ — 59° region since it is charaterized by an inverse P- 
Cygni profile overlaid on a broad component (see Section 4.2.5). 

N 2 H + and NH 3 . While the N 2 H + line-profile looks quite 
Gaussian, and thus used to set the source velocity, NH 3 shows an 
asymmetric line-profile which exhibits self-absorption (less evi¬ 
dent in 159-442), though not at the exact same velocity as absorp¬ 
tion in the water line (see Table 3). Absorption is also present in 
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mentioned, in the sources 130-304 and 130-327 the Gaussian 
component in absorption appears to be almost saturated. No 
line wings are detected above the noise level. In several sources 
(e.g., 130-304 and 130-327) the presence of multiple water 
absorption lines is evidently tracing the diffuse interstellar gas 
(which will be the subject of a separate paper, Persson et al., 
in prep.), that are likely a consequence of the larger amount of 
diffuse gas encountered along the £ = 30° direction. In Fig. 4 
we have already shown that the water line luminosity toward the 
£ = 30° sources is lower compared to that of the £ = 59° region. 

N 2 H + and NH 3 . As for the £ = 59° region, the N 2 H + line- 
profile looks Gaussian, and (when detected) it is also used to set 
the source velocity. In most cases the line-profiles of water and 
ammonia are consistent, i.e., they are either both in emission or 
absorption. Toward source 130-304 the NH 3 spectrum exhibits 
an inverse P-Cygni profile and the component in absorption 
looks almost saturated as in the water line. Toward source 
130-327 (and possibly 130-42), ammonia also shows absorption 
lines along the line of sight, as already observed in the water 
spectra. The maximum velocity difference between NH 3 and 
N 2 H + , when both lines are seen in emission, is - 0.9 km s -1 in 
source 130-42. 

We have compared the line widths of the NH 3 and N 2 H + 
lines, as shown in Fig. 5. Except for source 130-304, which 
shows an inverse P-Cygni profile in the ammonia spectrum, in 
all other sources we note that AV(NH 3 )> AF(N 2 H + ) including 
all four detected sources in the £ = 59° region. We also note that 
AVTNH i) tends to be somewhat higher in the £ — 59° sources. 
Since molecular outflows are the most distinctive feature of the 
detected £ = 59° sources, this is an indication that the line 
width of NH 3 is being broadened by these outflows. Likewise, 
the N 2 H + emission must then originate in a volume where the 
molecular gas is not significantly affected by the molecular out¬ 
flows. 


Fig. 5: Top. Line width of the NH 3 line vs. line width of the N 2 H + 
line, when both are detected in emission. The dashed line repre¬ 
sents the AT(NH 3 ) = AV4N4H' ) locus. Bottom. Line width of 
the broad component of the water line profile vs. the line width 
of the medium component for the £ = 59° sources. Symbols are 
as in Fig. 2. 


the NH 3 spectrum of 159-444. Finally, we note that one of the 
Gaussian components of the NH 3 line toward sources 159-441 
and 159-442 has a line width of 10.1 km s~ ] and 6.4 km s ' 1 (see 
Table 3). This suggests that the outflows observed in the water 
line are also partially detected in the ammonia line, at least in 
source 159-441. 

4.2.2. Line profile components in the £ = 30° region and 
comparison of line widths 

The nature of the observed lines toward this region is sig¬ 
nificantly different from that of the sources in the £ = 59° 
region. The line intensities are in general lower compared to 
those measured toward the £ = 59° region, which is likely 
a consequence of the average larger distance to the £ = 30° 
sources. 

H 2 0. The water line is detected both in emission and ab¬ 
sorption, but the line-profiles are much more varied. As already 


4.2.3. Narrow component 

The general appearance of the narrow component, mostly in ab¬ 
sorption, indicates a high H 2 0 column density and a low excita¬ 
tion temperature, suggesting that the absorption is being caused 
by the outer envelope and ambient cloud. Source 130-376 is the 
only one showing a clear narrow component in emission, which 
could be related to temperature, as appearance in emission re¬ 
quires an excitation temperature above the continuum level. We 
note, in fact, that 130-376 has the lowest continuum tempera¬ 
ture (23 mK, see Table 2) among sources in the £ = 30° region 
with a detectable continuum emission. As mentioned earlier, in 
the other £ = 30° sources the water line, either in absorption or 
emission, is broader than - 5 km s _1 and has thus been classified 
as a medium component. 

As already noted by Kristensen et al. (2012), because the ab¬ 
sorption of the narrow component is seen against both the out¬ 
flow and continuum emission, the absorbing layer must be lo¬ 
cated in front of both the emitting layers, i.e., the outflows must 
also be embedded. Figure 6 shows that the narrow component 
is either centred at the source velocity (offset 4 1 km s _1 ) or 
blueshifted (130-512) relative to the systemic velocity, as given 
by N 2 H + . In most £ = 59° sources, however, the narrow com¬ 
ponent is also slightly blushifted relative to N 2 H + emission, 
suggesting that the envelope, where the absorbed narrow com¬ 
ponent originates, is expanding. In addition. Fig. 6 shows that 
the narrow component is also blushifted relative to both the 
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Fig. 6: Central velocities of the narrow (V nar ) and medium (V me d) 
Gaussian components of the observed FFO lines compared with 
each other and with the systemic velocity from N 2 HL Symbols 
are as in Fig. 2. 


medium and broad (not shown) components. The only source 
where FnarlFLO) is red-shifted compared to both FmedfFLO) and 
Fbroadd^O) ' s 159-444, which has an inverse P-Cygni profile, 
which is instead typical of an infalling envelope. 

Contrary to van der Tak et al. (2013) we do not find any cor¬ 
relation between the source mass (which is mostly due to the 
envelope contribution, being derived from the Hi-GAL observa¬ 
tions of cold dust emission) and the line width of the narrow 
component. In Section 3.4 using RADEX we have estimated the 
0 -H 2 O column density of the envelope gas responsible for the 


observed absorption toward all £ = 59° sources and toward the 
£ = 30° sources 130-304,130-327 and 130-512 (see also Table 3). 


4.2.4. Outflows and shocked gas 

The H 2 O line profiles are characterized by two additional and 
broader components: a broad emission component, always in 
emission and visible only towards the £ = 59° sources, and a 
medium absorption or emission component seen in almost all 
sources. The broad and medium line components have widths 
exceeding the line width typical of the more quiescent envelope 
(e.g., Jprgensen et al. 2002) and are thus likely related to a proto- 
stellar outflow. 

In the £ = 59° sources, both Gaussian components are re¬ 
quired to fit the line profile. The medium component has al¬ 
ready been observed by Kristensen et al. (2012), toward a sam¬ 
ple of low-mass proto-stellar objects, by van der Tak et al. 
(2013) toward a sample of high-mass star forming regions, and 
also by Herpin et al. (2012) toward the massive protostar W43- 
MM1. Figure 6 shows that the medium component in most 
cases is centred at the source velocity (offset ;$ 1 kms" 1 ), with 
VjnedfFLO) - Vjsi-tNLH 1 ) having positive values in almost all 
cases. We already noted in the previous section that the nar¬ 
row component is blueshifted relative to the medium component, 
van der Tak et al. (2013) found that their medium component 
(which they call “narrow outflow component”) usually appears 
in absorption and is instead blueshifted relative to the narrow 
component (which they call “envelope component”). However, 
it should be noted that their decomposition of the line profiles is 
based on higher-excitation water lines. 

Both Kristensen et al. (2012) and van der Tak et al. (2013) as¬ 
sociate the origin of the medium component to the region, closer 
to the central protostar, where the gas is currently being shocked, 
whereas they associate the broader component with the swept- 
up gas of the molecular outflow. However, Herpin et al. (2012) 
report that their medium velocity component is likely due to a 
combination of turbulence and infall. A more detailed analysis 
of the different velocity components in shocked molecular gas 
has been carried out, for example, by Busquet et al. (2014) in 
the L1157-B1 shock region. They performed a multi-line H 2 O 
survey probing a wide range of excitation conditions. In addi¬ 
tion to the component in absorption (probably due to the pho¬ 
todesorption of icy mantles at the edge of the cloud), two gas 
components, connected with a J-type shock associated with a 
jet and with different excitation conditions, were found. The 
first component is a relatively extended (~ 10" corresponding to 
~ 2 x 10 4 AU), warm (~ 200 K), and dense (~ 10 6 cm -3 ) region, 
and the second is a more compact (less than 5" or ~ 10 4 AU), 
hot (~ 1000 K) and tenuous (~ 10 3 - 10 4 cm -3 ) component. The 
water abundance is very high in the hot gas (~ 10 4 ) and it de¬ 
creases down to ~ 10 6 in the warm component. Both warm and 
hot components have their peak at low-velocities, with the hot 
region dominating the total flux at the high end of the velocity 
range. 

In our work, and with only the data presented here, it is not 
possible to infer the exact physical origin of the medium compo¬ 
nent of the H20(1 io - loi) transition. However, the correlation 
between the line widths of the broad and medium components 
in the t = 59° sources, shown in the bottom panel of Fig. 5, 
suggests an outflow or shock origin for the medium component 
as well. We also note that Kristensen et al. (2012) speculate 
that the phenomenon producing the masers should also have a 
thermal component associated with it, which then produces the 
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medium component. We analyze the association with masers in 
Section 4.3. 

As far as sources in the £ = 30° region are concerned, no 
broad component has been detected. Observationally, this may 
be a consequence of the larger distances of these sources com¬ 
pared to the £ = 59° region. Alternatively, the viewing angle 
of the molecular outflows and/or shocks may also determine 
whether a source shows a broad, or both a broad and medium 
component. Because of the small sample size of £ = 30° sources, 
such an effect could be statistically significant, although in the 
£ = 59° region, with a comparable sample size, all sources show 
a broad component. An alternative, more physical reason for the 
lack of a broad component toward the £ = 30° region could be a 
different evolutionary phase (see Section 4.4). 

As far as NH 3 and N 2 H + are concerned, the line width of 
the ammonia line is likely affected by the presence of molecular 
outflows. We note, in particular, that source 159-441 is the only 
one where the NH 3 profile shows the same three components as 
the water line. While the water narrow component is generally 
blueshifted relative to the medium and broad components, all 
ammonia components have remarkably similar velocities. Three 
more sources show a NH 3 narrow and medium component (or 
two narrow components): 159-442, 159-444, and 130-304. The 
two ammonia components have very similar velocities also in 
159-442. In the other two sources, 159-444, and 130-304, how¬ 
ever, inverse P-Cygni profiles are observed and will be discussed 
in the next Section. 

4.2.5. Regular and inverse P-Cygni profiles 

In a total of three cases, we have observed a regular or inverse 
P-Cygni profile, which are usually associated to envelope expan¬ 
sion and infall, respectively. An inverse P-Cygni profile is char¬ 
acterised by an absorption feature which is red-shifted relative 
to an emission feature, and is typically caused by a velocity and 
excitation gradient in the gas falling in toward a central object. 
Two sources have inverse P-Cygni profiles, 130-304 and 159-444. 
In the case of source 159-444, the inverse P-Cygni profile is visi¬ 
ble in the spectrum of water as superposed on the broad outflow 
component. In this source a red-shifted absorption is also visible 
in the ammonia spectrum. Toward source 130-304 the inverse P- 
Cygni profile is only detected in the NH 3 spectrum, whereas the 
water spectrum shows multiple absorption lines from the diffuse 
interstellar gas, and at the source systemic velocity the water line 
profile appears almost completely self-absorbed (Section 3.4.2). 
A recent example of multiple ammonia lines, observed with 
Herschel toward the massive star forming region G34.3+0.15, 
showing inverse P-Cygni profiles can be found in Hajigholi et al. 
(submitted). The only regular P-Cygni profile we have detected, 
toward source 130-512, is also superposed on what looks like an 
outflow/shock component. These three profiles are thus all dif¬ 
ferent at some level and suggest a great variety of kinematic and 
excitation conditions. 

In order to give the simplest possible description of a con¬ 
tracting or expanding system that can reproduce the observed 
line profiles, Di Francesco et al. (2001) used a simple model con¬ 
sisting of two slabs moving towards each other with a continuum 
source between them. This model was a modification of a pre¬ 
vious “two-layer” code used by Myers et al. (1996). This model 
is appealing because of its simplicity, but on the other hand it 
gives only the crudest picture of the optical depth and excitation 
temperature of the absorbing and emitting gas. In addition, the 
model uses a relatively large number of free parameters (eight 
in total, if none is held fixed), including the excitation temper- 
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Fig. 7: From top to bottom, inverse P-Cygni profiles detected 
toward sources 159-444 (FFO) and 130-304 (NH 3 ), and regular 
P-Cygni profile detected toward 130-512 (HiO). The infall (or 
expansion, in the case of 130-512) model is shown on top (red, 
dashed line). The broad outflow component has been subtracted 
from the spectrum of source 159-444, whereas a narrow absorp¬ 
tion component has been added back to the spectrum of 130-512. 


ature of the front/rear layers and the continuum source (T/, T r , 
T c ) the filling factor of the absorbing source with respect to the 
telescope beam (®), the peak optical depth of the front and rear 
layers (t c ), the turbulent velocity dispersion in each layer (cr), 
the infall velocity (Vj n ) and the systemic velocity (Vj sr ). 
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Table 6: Best-fit parameters of the infall/outflow modelling. See text for a description of the parameters used in the fit. 


Source 

Species 

Fin 

[km s -1 ] 

cr 

[km s -1 ] 

Vi*/ 

[km s~'] 

To 

<h 

r j a 

1 c 

[K] 

Tf a 

[K] 

T r 

[K] 

159-444 

h 2 o 

0.5 

1.5 

34.3 

1.7 

0.14 

37.0 

4.0 

17.2 

130-304 

nh 3 

2.5 

1.7 

90.7 

3.1 

0.01 

350.0 

4.0 

10.5 

130-512* 

h 2 o 

-2.1 

1.6 

61.4 

1.3 

0.01 

63.0 

4.0 

11.0" 


Notes. (a> Fixed parameter, see text . ,b) Tentative result. 


Therefore, applying the simple two-slab model described 
above to the two sources with inverse P-Cygni profiles, we ob¬ 
tain the best-fit infall models shown in Fig. 7, with the corre¬ 
sponding best-fit parameters listed in Table 6. In the case of 
source 159-444, the broad outflow component (and an additional 
narrow absorption feature, blue-shifted with respect to the in¬ 
verse P-Cygni profile, see Table 3) have been subtracted to the 
original spectrum before applying the model. For both sources 
the value of T c has been fixed to that shown in Table 2. We 
have also fixed the value of Vi sr to the velocity of NAH 1 (see 
Table 3). In section 4.2.3 we suggested that the gas front layer 
may have a low excitation temperature and thus we fixed it to 
the approximate intermediate value {Tf =* 4 K) shown by Myers 
et al. (1996) in their Fig. 1. As a comparison Hajigholi et al. (sub¬ 
mitted) find an excitation temperature of ^ 7 - 8 K in the absorb¬ 
ing envelope of source G34.3+0.15. 

By comparing different infall model profiles to the obser¬ 
vations we find that Vj n and cr can vary by as much as 0.2 - 
0.5 kms -1 , and are not very sensitive to Tf. However, these tests 
do not take into account the uncertainty on the beam-filling fac¬ 
tor, used to convert the antenna temperature to line brightness 
temperature (see Section 3.4), which may also affect the best-fit 
results. Thus, our results show that while the infall velocity to¬ 
ward 159-444 is comparable to the values found by Di Francesco 
et al. (2001) and Kristensen et al. (2012) (who analyze low-mass 
protostars), the values of cr (and the value of V 7 ;,, toward 130-304) 
are somewhat higher, since these authors find Vj n , <x < 1 kms -1 . 
Our estimated values of cr are consistent with those found by 
Hajigholi et al. (submitted) toward G34.3+0.15. Some of the 
high-mass protostars studied by van der Tak et al. (2013) also 
show regular and inverse P-Cygni profiles, but they do not an¬ 
alyze the line profiles and thus we cannot directly compare our 
results. Instead, we can compare our results with those of Herpin 
et al. (2012) toward the massive protostar W43-MM1, where 
they derive infall velocities as high as 2.9 km s 1 and turbulent 
velocities ^ 2 km s 1 . These values are much much more similar 
to our estimates of cr and of the infall velocity toward 130-304. 

As far as the only source, 130-512, with a regular P-Cygni 
profile in our sample, this source shows a narrow absorption 
component at a velocity of 64 km s 1 . Therefore, assuming that 
this additional absorption feature is not related to the overall 
P-Cygni profile, and in order to reconstruct the hypothetical P- 
Cygni profile we added back this component to the water spec¬ 
trum. The two-layers model described earlier can also be applied 
to the case of expanding gas, the only change in the model is that 
the infall velocity becomes negative. We show our best-fit result 
in the bottom panel of Fig. 7, obtained by fixing the value of T r , 
which would otherwise diverge. Clearly, the fit is not very good 
and it is difficult to determine whether this may be due to con¬ 
tamination of the line profile, up to the level of possibly resulting 


in a fake P-Cygni profile. For this reason our results in Table 6 
are labelled as tentative. 

While the infall and turbulent velocities are distance- 
independent and, at least in this toy model, do not depend on the 
mass distribution within the source, determining the mass infall 
rate would require the knowledge of the characteristic distance 
of the infall zone or the mass undergoing infall. If the envelopes 
are in free fall, the velocity at a radius R m in the envelope is given 
by Vin = y/2GM enc /R{ n where M enc is the mass enclosed within 
the radius R m . Assuming spherical infall symmetry, the mass in¬ 
fall rate can be estimated as (e.g., Kristensen et al. 2012): 

M = 4nfim H n m (2GM cnc ) 2 V~ 3 (4) 

with n m the density at the infall radius, mu and p = 2.8 the 
hydrogen atom mass and mean molecular weight, respectively. 
Therefore, in order to estimate M it is necessary to be able 
to determine either M enc or the infall radius R m . We do not 
have a reliable distance determination for source 159-444 (see 
Table 2), and thus we cannot estimate any mass. In the case 
of source 130-304 we do have a distance estimate and a mass, 
but the large mass and relatively large distance (see Table 2) 
suggest that the estimated mass represents a much larger mass 
than that enclosed within the infalling region. As an example, if 
we used a typical infall radius estimated by Herpin et al. (2012) 
(R m ~ 6 x 10 17 cm) then in source 130-304 we would get a mass 
infall rate M ~ 5 x 10 2 M 0 yr ] , using as n m the density given in 
Table 2 and V m from Table 6. Not surprisingly, M is comparable 
to the values estimated by Herpin et al. (2012) toward the mas¬ 
sive proto-stellar clump W43-MM1, and is much higher than the 
infall rates determined by Kristensen et al. (2012) toward low- 
mass proto-stellar sources. 

4.3. Association with maser emission 

We searched for masers associated with the detected sources 
(within a 30" radius) and found a total of six sources associated 
with maser emission: all four detected sources in the ( = 59° 
region and two more in the £ = 30° field. Sources 159-441 and 
159-445, being very close, are actually associated with the same 
methanol maser. All sources, except for 159-442, are associ¬ 
ated with a methanol 6.7 GHz maser. Source 159-442 is the only 
one associated with a 22 GHz water maser (B 193852.5+235736, 
Palla et al. 1991), while source 130-304 is also associated with 
an OH maser. 

The methanol masers associated with sources 159-441/159- 
445 and 159-444 were found by comparing our observations with 
those carried out by Olmi et al. (2014b) toward the ( - 30° and 
{ = 59° regions in search of 6.7 GHz CH 3 OH and 6.0 GHz OH 
masers. The masers associated with 159-441/159-445 and 159- 
444 are G59.63-0.19 and G59.78+0.63, respectively, and they 
have peak flux densities of 0.58 and 0.03 Jy. As discussed by 
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Olmi et al. (2014b), these two methanol masers belong to a new 
class of low-brightness masers that could represent an earlier 
stage of evolution. 

The methanol masers have velocities which agree well with 
the velocities observed in our detected sources (differences are 
$1-2 km s _1 ). However, the velocity of the 22 GHz water maser 
(Vi sr 50kms _1 ) associated with source 159-442 significantly 
differs from the central velocity of any of the Gaussian com¬ 
ponents of its HjO line profile (Vi sr ^ 35 - 37kms~ 1 ). If we 
assume that the physical association between these two sources 
is indeed real, given the small angular separation between them 
(~ 4 arcsec), then one might speculate that the water maser is 
originated somewhere along the path of the outflow. In fact, the 
broad Gaussian component in 159-442 has a line width that en¬ 
compasses the velocity of the water maser. 

4.4. Evolutionary trends 

The relative intensities of the various components used to de¬ 
compose the line profiles of the water spectra have already been 
used in an attempt to identify possible evolutionary sequences. 
For example, Kristensen et al. (2012) associate the disappear¬ 
ance of the broad component with later evolutionary stages in 
their low-mass protostars sample, van der Tak et al. (2013) then 
propose that the appearance of the narrow component (or “en¬ 
velope” component) in emission is an indicator of a warmer en¬ 
velope, and thus associated to later evolutionary stages. We note 
that while in all £ — 59° sources the narrow component appears 
in absorption, toward the £ — 30° region sources 130-43 and 
130-376 are characterized by a narrow component in emission 
(see Table 3; see also Section 3.2 regarding the water detection 
toward 130-43). In source 130-512 the narrow component is in 
absorption, but this is a more complex source with a regular P- 
Cygni profile, as discussed in Section 4.2.5. 

An alternative method to analyze the evolutionary stages 
of our sources makes use of the evolutionary tracks previously 
shown in Fig. 2. These tracks were determined following a 
method similar to that discussed by Molinari et al. (2008), who 
proposed an evolutionary sequence for high-mass protostars in 
terms of two parameters: the envelope mass and the bolomet- 
ric luminosity. These authors suggest an evolutionary sequence 
which is concentrated in two main phases: protostars first accrete 
mass from their envelopes, and later disperse their envelopes 
by winds and/or radiation. In Fig. 2 the calculated evolutionary 
tracks thus rise upward almost vertically in luminosity during 
the accretion phase, and then proceed horizontally to the left (to 
lower masses) during the envelope dispersal phase. The posi¬ 
tions of our sources in this plot suggest that some (but not all) of 
the sources with at least one line detection are, for similar ranges 
of mass, in a more advanced stage of evolution. This is certainly 
the case for the detected £ = 59° sources with a well defined dis¬ 
tance (159-441 and 159-445), while for sources detected toward 
the £ — 30° region a preferential evolutionary phase is not appar¬ 
ent. None of these sources, however, seem to belong to the late 
phase of envelope dispersal. 

There are other indicators that can be used to analyze the 
relative evolutionary phase of the £ = 30° and £ = 59° re¬ 
gions. For example, from the discussion about T c in Section 4.1 
and Table 5, there is a trend 10 of higher T c toward the de¬ 
tected £ = 30° sources, suggesting they may be in a more ad¬ 
vanced stage of proto-stellar evolution, compared to sources in 


10 We use the word “trend” because of the large scatter shown by the 
T c values in Table 5. 


the £ = 59° region. Furthermore, as already noted in Section 3.3, 
the typical water line luminosity toward the £ = 30° sources is 
lower compared to that of the £ = 59° region. Kristensen et al. 
( 2012 ) find a similar decrease in water emission between their 
Class 0 and Class I low-mass protostars. Finally, as stated in 
Section 4.3 two of the £ — 59° sources are associated with low¬ 
brightness methanol masers that could indicate an earlier phase 
of evolution. Source 130-304 is the only one associated with an 
OH maser, which is sometimes described as characterizing later 
evolutionary stages (Breen et al. 2010) compared to 6.7 GHz 
methanol masers. Therefore, a firm conclusion about the spe¬ 
cific evolutionary phase of the observed sample is limited by the 
small number of sources. However, various indicators tentatively 
suggest that the detected sources toward the £ = 30° region are 
in a later evolutionary phase compared to those in the £ — 59° 
region. 

5. Summary and conclusions 

Using the HIFI instrument of Herschel we have observed the 
0 -NH 3 G 0 - 0 0 ) (572 GHz), o-H 2 0(1io - lot) (557 GHz) and 
N 2 H + (6 - 5) (559 GHz) lines toward a sample of 52 high-mass 
starless and proto-stellar clumps selected from the Hi-GAL sur¬ 
vey. The target sources are concentrated in two galactic fields, 
at £ = 30° and £ = 59°, and were observed during the Science 
Demonstration Phase of Herschel. At least one of the three se¬ 
lected molecular lines was detected in 4 out of 35 sources toward 
the £ = 59° region, and 7 out of 17 objects toward the £ = 30° 
region. 

Almost all water spectra can be described as the sum of a nar¬ 
row (FWHM- 3-5kms -1 ), a medium (FWHM- 6 -lOkms- 1 ) 
and a broad (up to FWHM- 30km s 1 ) velocity component. All 
detected sources of the £ — 59° region have a strikingly similar 
water line-profile, characterized by the broad Gaussian compo¬ 
nent associated with the presence of molecular outflows, which 
are also observed in the NH 3 line, at least toward one source. 
Therefore, in these sources there is strong evidence that water 
line emission is linked to the presence of shocked molecular gas. 

The situation is very different in the £ = 30° region, with 
no two profiles being identical, and with no clear detection of a 
broad component in the water spectra. We find that the contin¬ 
uum temperature, T c , is correlated with the source bolometric lu¬ 
minosity. We also find an indication that the water line luminos¬ 
ity increases with T c . The detected sources in the £ = 59° region 
have both a higher water line luminosity and water abundance, 
compared to sources detected in the £ = 30° region. Although 
we have few sources, this trend constitutes additional support to 
a scenario where the water abundance is linked to shocks. 

In terms of the gas kinematics, the general appearance in 
absorption of the narrow Gaussian component indicates a high 
H 2 0 column density and a low excitation temperature, suggest¬ 
ing that the absorption is being caused by the outer envelope and 
ambient cloud. The correlation between the line widths of the 
broad and medium components in the £ - 59° sources suggests 
an outflow or shock origin for the medium component. In addi¬ 
tion, we also detect in a few sources inverse and regular P-Cygni 
profiles, though we do not find any link with other source param¬ 
eters. A firm conclusion about the specific evolutionary phase of 
the observed sample is limited by the small number of sources. 
However, various indicators tentatively suggest that the detected 
sources toward the £ = 30° region are in a later evolutionary 
phase compared to those in the £ — 59° region. 
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Appendix A: Spectra of sources in the t = 30° 
region 

All spectra are in T* A units and have been smoothed to a spectral 
resolution of 1 kms -1 , unless noted otherwise. 





























SSB Antenna temperature T. [K] SSB Antenna temperature 
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Appendix B: Spectra of sources in the £ = 59° 
region 

All spectra are in T* A units and have been smoothed to a spectral 
resolution of 1 kms -1 , unless noted otherwise. 







































SSB Antenna temperature T. [K] SSB Antenna temperature 



























































































































